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STUDIES  OF  THE  ELECTROMAGNETIC  CHARACTERISTICS 


OF  MOVING  IONIZED  GASES 


The  Project  Scientist  for  this  task  is  Dr.  Manlio  Abele. 

The  investigations  of  the  microscopic  and  macroscopic 
properties  in  a  weakly  ionized  gas  reported  in  Chapters  II  and  III 
of  this  summary  report  were  carried  out  under  the  direction  of 
Professors  Piero  Caldirola  and  Luigi  NapoLitano,  respectively, 
and  under  the  general  direction  of  Prof.  Piero  Caldirola. 

The  studies  of  the  propagation  of  electromagnetic  waves 
in  the  shock  tube,  reported  in  Chapter  IV,  were  carried  out  by 
Drs.  Frank  Lane  and  Gino  Moretti. 

The  experimental  work  reported  in  Chapter  V  has  been 
carried  out  under  the  direct  supervision  of  Dr.  Abele  assisted  by 
Messrs.  Roger  Tomboulian  and  Myron  Wecker. 


CHAPTER  I  ' 


INTRODUCTION 

This  report  contains  preliminary  information  on  the  theoretical  analysis 
of  the  electromagnetic  properties  of  a  non-uniform  ionized  gas  and  the  des¬ 
cription  of  the  facility  to  be  used  for  the  experimental  investigation  of  the 
propagation  of  an  electromagnetic  wave  in  the  non-uniform  gas. 

It  is  well  known  that  in  a  first  order  approximation,  the  propagation 
of  an  electromagnetic  wave  may  be  analyzed  with  a  simplified  scheme  in 
which  the  local  properties  of  the  ionized  gas  are  given  by  the  plasma  frequency 
and  the  constant  collision  frequency.  Even  in  this  case  the  solution  of  two- 
or  three-dimensional  propagation  problems  may  present  tremendous  mathe¬ 
matical  difficulties.  On  the  other  hand,  in  a  general  case  for  a  non-uniform 
gas,  the  effort  required  to  solve  these  problems  is  not  justified  due  to  the 
oversimplified  scheme  of  the  electrical  properties  of  the  gas.  When  the 
change  of  the  macroscopic  thermodynamic  properties  of  the  gas  becomes 
important  in  a  length  of  the  same  order  of  magnitude  of  a  wa\  e  length  of  the 
electromagnetic  field,  it  becomes  netessary  to  analyze  the  possible  effects 
of  additional  phenomena  such  as  the  transport  process,  the  anisotropic  be¬ 
havior  of  the  electrical  properties,  the  electrical  polarization  of  the  medium, 
etc*  Furthermore,  the  coupling  between  transverse  oscillations  and  coherent 
longitudinal  e  1  e  c  t  r  o  n  os c illations  may  become  a  significant  factor  leading 
to  a  substantial  increase  of  the  dissipation  of  electromagnetic  energy  in 
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thermal  energy  in  the  ionized  gas.  If  the  relative  importance  of  these  phen¬ 
omena  could  be  properly  evaluated  from  the  theoretical  standpoint,  then  it 
would  be  possible  to  define  the  set  of  parameters  which  govern  the  behavior 
of  an  ionized  gas  for  each  particular  case.  The  major  difficulty  of  a  theo¬ 
retical  approach  is  due  to  the  fact  that  it  requires  detailed  information 
about  the  microscopic  processes,  which  at  the  present  time  is  rather  in¬ 
complete,  particularly  for  the  case  of  weakly  ionized  air  as  encountered  in 
fl  ight  conditions. 

Thus  a  systematic  experimental  investigation  becomes  of  primary  im¬ 
portance  in  order  to  obtain  direct  information  about  the  electromagnetic 
properties  and  it  must  be  complemented  by  a  theoretical  analysis  of  the 
microscopic  processes  which  occur  in  a  non-uniform  gas.  Furthermore, 
a  theoretical  research  has  to  be  developed  in  order  to  establish  the  correla¬ 
tions  between  the  experimental  data  and  the  results  of  the  theoretical  analysis 
of  the  physical  properties  of  the  gas. 

The  experimental  investigation  of  the  local  electromagnetic  properties 
of  a  non-uniform  gas  by  means  of  the  study  of  the  propagation  of  an  electro¬ 
magnetic  wave  involves  a  huge  amount  of  difficulties.  An  overall  measurement 
of  the  total  phase  shift  and  attenuation  of  the  wave  across  the  gas  in  a  given 
direction  would  not  provide  enough  information  about  the  local  properties. 

On  the  other  hand,  a  direct  measurement  of  the  distribution  of  electromag¬ 
netic  field  must  be  obtained  without  perturbing  both  the  flow  field  and  the 
electromagnetic  wave.  Consequently,  the  geometry  of  the  channel  where 
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the  propagation  is  analyzed,  must  be  such  that  the  necessary  information 
may  be  obtained  from  a  set  of  measurements  performed  at  the  boundary  of 
the  channel  itself.  In  the  case  of  a  one-dimensional  non-uniform  flow,  these 
considerations  suggest  the  advantage  of  using  the  channel  containing  the 
ionized  gas  as  an  element  of  wave  guide  where  certain  basic  conditions  about 
the  distribution  of  the  electromagnetic  field  may  be  predicted  provided  that 
the  proper  range  of  frequencies  is  selected  for  the  electromagnetic  field. 

By  using  a  shock  tube  to  produce  the  ionized  gas,  the  most  simple  experi¬ 
mental  arrangement  corresponds  to  the  propagation  of  the  electromagnetic 
wave  in  a  section  of  the  driven  paft  of  the  shock  tube  itself.  A  non-uniform 
distribution  of  the  electrical  properties  of  the  gas  may  be  achieved  by  es¬ 
tablishing  a  temperature  gradient  in  the  driven  section,  as  discussed  in 
detail  in  Chapter  V  of  this  report.  If  the  flow  were  uniform  in  a  cross-section 
of  the  wave  guide,  with  a  suitable  distribution  of  microwave  detectors  located 
at  the  wall  of  the  wave  guide,  we  would  obtain  complete  information  about 
the  properties  of  propagation  of  the  electromagnetic  field.  In  turn,  this  in¬ 
formation  coupled  with  the  theory  of  propagation  in  a  wave  guide  leads  to  the 
correlation  of  the  measured  quantities  with  the  theoretical  values  of  the  phy¬ 
sical  properties  of  the  gas.  Of  particular  interest  is  the  correlation  between 
experimental  and  theoretical  values  when  the  gradient  of  electron  density  is 
chosen  in  such  a  way  that  the  electromagnetic  wave  cannot  propagate  beyond 
a  given  cross-section  of  the  ionized  gas.  The  position  of  this  critical  cross- 
section  of  the  ionized  gas  and  the  Amount  of  reflection  of  the  electromagnetic 
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wave  at  this  section  should  be  very  sensitive  to  a  small  change  in  the  electrical 
parameters  of  the  gas.  Thus  we  may  expect  the  maximum  detectable  effect 
of  a  departure  of  the  actual  physical  properties  from  the  theoretical  values. 
This  critical  cross-section  moves  along  the  shock  tube  with  the  velocity  of 
the  ionized  gas.  A  microwave  detector  located  in  the  region  of  propagation 
before  the  arrival  of  the  critical  cross-section  behaves  like  the  detector  of 
a  standard  standing  wave  ratio  meter.  After  the  instant  of  time  at  which  the 
critical  cross-section  reaches  the  position  of  the  detector,  the  detector  will 

indicate  the  damping  of  the  electromagnetic  field  in  the  region  of  total  atten¬ 
uation. 

The  position  of  the  initial  cross-section  may  be  controlled  with  an 
externally  applied  constant  magnetic  field.  With  an  axial  distribution  of  mag¬ 
netic  field  which  should  not  perturb  the  gross  properties  of  the  flow  field,  we 
introduce  an  additional  parameter  that  may  be  extremely  helpful  in  obtaining 
the  information  about  the  new  properties  of  propagation  which  are  created  in 
the  wave  guide.  If  the  mi  croscopic  properties  of  the. collision  processes  in 
the  gas  are  not  strongly  affected  by  the  external  magnetic  field,  the  anisotropic 
electrical  properties  of  the  ionized  gas  may  be  predicted  on  the  basis  of  the 
solution  of  the  equations  which  govern  the  macros copic  motion  of  the  electrons. 

A  careful  analysis  of  this  particular  problem  from  both  the  microscopic  and 
macroscopic  point  of  view  is  one  of  the  objectives  of  the  theoretical  investi¬ 
gation  outlined  in  this  report.  ' 
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The  assumption  of  a  one -dimensional  flow  is  verified  in  the  shock  tube 
only  in  a  first  order  approximation.  Actually,  the  boundary  layer  behind  the 
shock  wave  introduces  a  non-uniform  distribution  of  the  electrical  properties 
in  a  cross-section  of  the  shock  tube.  The  propagation  jncde  chosen  for  the 
wave  guide  is  the  fundamental  T.  E.  mode,  and  the  frequency  of  the  electro¬ 
magnetic  wave  is  such  that  this  is  the  only  mode  that  can  be  propagated. 

For  these  conditions  the  thickness  of  the  boundary  layer  where  a  strong  change 
of  the  electrical  properties  may  be  expected,  is  a  small  fraction  of  the  wave 
length.  Consequently,  the  non-uniform  transverse!  distribution  of  the  flow 
field  will  appear  as  an  average  effect  on  the  axial  propagation  of  the  electro¬ 
magnetic  wave.  This  effect  must  be  analyzed  in  order  to  establish  to  what 
extent  it  may  be  neglected.  A  preliminary  calculation  is  presented  in  this 
report  with  a  simplified  scheme  which  should  give  the  order  of  magnitude 
of  the  effect  to  be  expected  in  the  experimental  conditions  of  the  test  section 
of  the  shock  tube.  It  is  necessary  to  point  out  that  the  boundary  layer  may 
introduce  additional  difficulties  corresponding  to  polarization  effects  and  dif¬ 
fusion  processes  due  to  the  strong  transverse  gradients,  and  a  detailed  in¬ 
vestigation  of  these  effects  may  become  necessary  during  the  development  of 
this  research. 

On  the  basis  of  the  aforementioned  considerations  and  the  particular 
experimental  technique  to  be  used  in  the  measurement  of  the  electrical  prop¬ 
erties  of  the  non-uniform  flow,  the  theoretical  research  program  is  divided 
in  two  main  parts.  One  part  is  devoted  to  the  analysis  of  the  microscopic 
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phenomena  and  the  calculations  of  the  parameters  which  define  the  physical 
macroscopic  properties  of  the  non-uniform  gas.  As  shown  in  Chapter  II, 
the  statistical  approach  assumes  a  weakly  ionized  gas  in  such  a  way  that  the 
distribution  function  for  the  molecular  species  is  independent  of  the  distribu¬ 
tion  function  for  the  electron  gas.  Thus  the  electric  properties  of  the  gas 
are  evaluated  in  terms  of  the  thermodynamic  properties  of  the  non-uniform 
flow  field.  This  assumption  is  certainly  valid  in  the  particular  condition  of 
the  gas  in  the  shock  tube. 

In  the  other  part,  the  theoretical  analysis  of  the  propagation  inside  the 
wave  guide  will  be  developed  with  different  model  distributions  of  the  macro¬ 
scopic  properties  of  the  ionized  gas  in  order  to  interpret  the  experimental 
results  and  to  establish  the  correlations  between  the  measured  electrical 

.V* 

properties  and  the  corresponding  theoretical  values. 
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CHAPTER  II 


ANALYSIS  OF  THE  MICROSCOPIC  PROPERTIES  IN  A  WEAKLY  IONIZED  GAS 

t.  INTRODUCTION 

The  main  object  of  the  work  presented  in  this  section  is  to  study  the  effects 
of  inhomogeneitie s  in  a  plasma  caused,  for  instance,  by  electric  fields,  magnetic 
fields  and  temperature  gradients. 

Let  us,  at  first,  consider  a  plasma  model  without  boundary  conditions,  ancj 
assume  that  the  electronic  density  is  much  smaller  than  the  molecular  density, 
(Lorentz  gas.)  In  doing  so,  we  may  overlook  in  our  model  the  presence  of 
the  ions.  Furthermore,  if  we  indicate  with  f(P,  v,  t)  and  F(P,  V,  t),  respec¬ 
tively,  the  electronic  and  molecular  distribution  functions,  we  may  overlook 
in  the  corresponding  system  of  two  coupled  Bolt/.inann  equations: 

Di  -  Jcv  (f)  +  Jen,(f) 

DK=  Jmni(F)4  J„U.(F> 

the  two  collision  operators  J  (with  respect  to  JC]K)  and  J  (with  respect  to 
Jmm).  A  list  of  symbols  is  given  at  the  end  ol  this  chapti  r 

Assuming  then,  the  following  system  for  the  temporal  buha\  i  >r  of  i  and  F, 

Df  =  Joal(f)  (1.1) 

DF=  Wn.  (I--) 

we  may  note  that: 

a.  From  the  equation  (l.J.)  the  temporal  behavior  of  the  molecules  is  entirely 
independent  of  the  presence  of  electrons,  while  the  electron  gas  depends  upon 
the  presence  of  molecules,  through  the  coltifeion  term  J, ,rrr» 


b.  The  temporal  behavior  of  the  electron  gas,  given  b\  equation  (l.l) 
is  affected  b^  the  molecular  gas,  as  we  shall  see  later,  only  through  the 
following  phenomenological  parameters:  molecular  density,  mean  velocity, 
and  mean  square  velocity. 

As  a  result,  it  is  not  necessary  in  solving  equation  (l.l)  to  integrate  the 
Boltzmann  equation  (I.Z): 

If  1-  -  Jmm( !• ) » 

because  we  only  need  the  knowledge  of  the  mentioned  phenomenological  parameters. 

“  •  ANALYSIS  of  THE  EQUATION  ( L  1_) 

Wo  shall  begin  by  expanding  tne  eh  ctronic  distribution  function  f(P,  v,t)  in 
spite  rical  harmoni  s,  writing: 


00 


f(p 


.  V.t)  )^ 


n  (P.’-.O  Cjinr  Pjrn(P,v.t) 


Jin 


(  c  ui<j  ) 


U.i) 


where; 

C jin=  -J  -os  n1<? 

sjm=  Vj  Sin  m°  ’ 

y  ^  (9)  being  the  associated  Legendre  functions: 

(9)  =  sin"1 0  P(m)  (cos  0). 

Jnl  J  ' 

We  may  write  the  expansion  (Z.i)  in  the  form: 

f  (v)  =  to(v)  +  4-  -  it  (v)  +X(i)  .  (Z.Z) 

where  f0  represents  the  isotropic  part  of  the  distribution  function  f,  fj  is  the 

coefficient  of  the  first  anisotropy  and  is  simply  the  remaining  part  of 

expansion  (Z .  1) . 
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Introducing  (1.  1)  into  the  Boltzmann  equation  (L.  1),  we  obtain  a  hierarchy 


of  equations,  which  is  easily  obtained  provided  we  know  the  result  of  the 
application  of  the  linear  collision  operator  J  over  each  term  of  the  expansion  (1 .  1). 

To  this  end,  we  may  note  that  the  operator  J,  in  an  approximation  for 
which  one  has  simply  m/M  =0  (perfect  Lorentz  gas),  possesses  the  following 
properties,  for  each  isotropic  function  a(>): 


J  (a) =  0 


J  Cjm>  =  ' 

VJ 

a  Cjm 

J  (.  Sjni)  =  - 

V 

a  S  • 

J 

jm 

In  the  imperfect  Lore 

ntz 

gas  scheme 

J  (a)  =  O^—) 
M 


we  shall  have: 


J  (Q  Cjm)  =  *  vj  a  Cjm  1  +  °z  ("KT* 


J(a  Sjm)  =  -  Vj  a  Sjm 


1+°3  <tt) 


(3.3) 


(3.4) 


where  we  have  used  the  Landau  symbol  O(x),  meaning  a  quantity  which  is 
simply  of  the  order  of  x. 

The  quantity  Oj  (— )  has  been  calculated  by  others  (Ref.  Z)  under  the 
M 

hypothesis  of  a  Maxwellian  molecular  distribution  F,  obtaining: 


J  (a)  =  —  _2_ 


l 


M 


v2  0v 


Vl 


(a  + 


kT 

mv 


Q  a 
dv 


(3.5) 


On  the  other  hand,  we  have  proceeded  to  the  evaluation  of  the  quantity 
0,(^)  keeping  an  entirely  arbitrary  (and  hence  generally  anisotropic)  molecular 
distribution.  The  result  of  our  calculation  is  given  by; 


T/  V  _  l  a 
J  (  a  )  -  — , - < 

v  dv  ' 


A  (v)  +  B(v)  v  •  V  +  C(v)  vv  :  VV  ^  , 


(-•  6) 


who  r  e: 


.  ,  >  _  m  i  ,  ,3  \  f)  n 

A  (\) - vj  v  a  +  v,v - —— 

M  j  d  v 


B(v) 

C(v) 


v  i  v'  a 


(3.7) 


,  l  t)  a  5  I  . 

(v,‘  '7  +  V-  L(vi 


1  .  1  e  V  2  2 

7  vz )  v  "  -  -  v 

L  o  0  V 


W  c  hdivi'  followed  Ci.  apmun  and  Cowling's  notation,  in  writing: 


vV  .  y  v 


where  V,  and  v-,  are  th<  first  two  relaxation  frequencies  (Ref.  1): 


TT 

r 

i 

1  sin  9 

l  -  P,  (COS  9)j 
»  ■* 

(-•  8) 


T(  r 

vz(\)  =  L  ttNv  J  sin  9  l  -  P2  (cos  9) 


(«.  v)  d9  , 


"(••Ml)  being  the  differential  cross  section  for  electron-molecule 

collision,  in  a  frame  of  reference  fixed  with  onr  of  the  two  particles. 

The  determination  of  J  (nCjm  ),  J  (  aSju  )  has  been  done  disre  garding 

the  02  ,  03  terms.  (We  will  give  a  bri  f  account  of  the  method  followed  in 

the  derivation  of  formula  (3.6)  later  in  this  chapter.)  Taking  into  account 

v  _  .  , 

that  fQ  is  an  isotropic  function  and  —  «  f  l  is  simply  a  linear  combination 

with  isotropic  coefficients  of  C  1Q,  Cn  ,  Su  ,  we  arc  now  able  to  calculate 

j(  f  \ .  J  (  .  f  )  ,  If  we  perform  a  furtlu  r  approximation*  forgetting  thi I'Xi'i ) 

\  O'  1  *  y  •—  1  1 

*  As  is  usual  in  this  kind  of  reasoning,  wc  may  evaluate  "a  posteriori"  the 
usefulness  of  this  approximation,  by  calculating  explicitly  the  difference  between 
the  results  obtained  and  those  based  on  the  successive  approximation,  or 
"a  priori"  we  may  introduce  some  plausible  arguments,  like  those  used  by 
Davydov  (Ref.  3)  and  Gurevitch  (Ref.  4). 


in  the  (Z.  J.)s  we  at  last  arrive  at  th<-  following 


D  f  =  J(i) 


=  [~2  ~  |a  *  Ci*v  :  vvj  -  v,  (-=-  .  Lj  .  (Z. 


9) 


where  A(v),  B(v)  and  C(  . )  are  given  by  (Z.7). 

Finally,  by  means  of  the  orthogonality  properties  of  the  spherical  harmonics 
functions  on  the  unity  sphere,  we  obtain  tin  following  system: 
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has  been  evaluated  by  noting  that  in  each  .  lectron-molecule  collision  the 
variation  of  the  modulus  v  of  electron  \eloe  ity  v  ,  is  given  by: 


S’-i  -  °|(t) 


lA 


(3.3) 


As  a  consequence,  the  integral  operator  (3.-1)  can  bi  put  in  a  differential 
form  following  the  same  technique  which  is  used  in  the  transformation  of  the 
Boltzmann  elastic  collision  operator  into  the  Fokker-Planck  form  (Refs.  5  andb). 


*  As  will  be  shown  later,  the  main  object  of  our  current  research  consists  in 
the  evaluation  of  stationary  or  asymptotic  solutions  of  the  system  (Z.  10). 
Simpler  systems,  in  which  we  may  give  evidence  to  particular  aspects  of 
inhomogeneities,  will  be  considered. 
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Let  uh  introduce  a  suitable  (but  largt  ly  arbitrary)  isotropic  function  ^  (v) , 


c cording  to  a  general  well  known  rule; 


Jj(  i)§(  )  <>Y  =  J 


§(v<)  -  §(v)j  gaF  a  («,  r)  di2  dV  dv  , 


(3.3) 


Furthermore,  retaining  only  terms  of  the  order  of  m/M,  we  have,  according 

to  (3 ,Z): 

r»X  i  •,  a-X 

(  3.4) 


Do  l  i  D2? 

$(V)  -0(v)  =  (V*-  v)  —  +  v)  — S- 

1  dV  t.  Bv- 


Inserting  (3.  1)  into  (3,. 4)  and  integrating  by  parts,  we  found: 


J(u)  =  —  ~ 
v“  dv 


>  ID  , 

-  v-  p  a  +  —  —  (  Vi) 

-  d  V 


(3.5) 


where 


(3.6) 


p(v>*  Jo  (8,g)  gF  (v'-  V)  di2  dV 

\(v  >  =Ju(4,g)  gF  (V-  v)2  d!2  dV  . 

Ii  a  way  consistent  with  the  approximations  used,  the  two  preceding  functions 
p  and  X  are  evaluated  at  the  first  order  in  the  quantity  m/M. 

The  rather  lengthy  calculations  are  omitted  here.  Assuming  a  quite 
arbitrary  molecular  distribution  function  F(V)  and  inserting  the  final  expressions 
for  p  and  X  into  (3.5),  we  obtain  the  formula  (a. 6).  We  note  explicitly  that, 
in  the  hypothesis  that  the  molecular  distribution  function  F  be  isotropic, 
formula  (Z.6)  reduces  to  (Z.5)  which  has  been  calculated  by  Chapman  and 
Cowling  (Ref.  Z) . 
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CONCLUDING  REMARKS 


A  preliminary  investigation  of  the  possibilities  of  solving  the  system 
(Z.  10)  is  now  in  progress, and  in  this  connection,  we  intend,  in  the  future, 
to  develop  the  following  points:  a.  To  first  integrate  the  system  (Z.  10)  for 
particular  cases  and  then,  if  possible,  in  general  form.  This  will  be  carried 
out  in  order  to  evaluate  the  electromagnetic  parameters  of  the  plasma  (taking 
inhomogeneities  into  account)  within  the  limits  of  the  approximation  used 
through  the  knowledge  of  the  electronic  distribution  function;  b  To  evaluate 
"a  posteriori"  the  usefulness  of  the  approximations  used;  and  c.  To  work  out 
some  possible  extensions  of  the  theory. 

a.  Being  aware  of  the  noteworthy  analytical  difficulties  of  our  system 
(Z.  10),  we  intend,  as  stated  before,  to  analyze  it  in  more  simple  cases,  at 
least  as  a  first  step.  For  instance,  if  we  take  E  =0,  H  =  0,  V  =0  and  we 
consider  a  stationary  stati  ,  the  system  (Z.  10)  reduces  to  one  equation  only 
in  the  isotropic  function  I0 ,  which  probably  can  be  integrated  in  closed  form. 
In  this  case,  with  the  further  hypothesis,  T  =  const.  ,  such  an  equation  is 
separable;  stating 

f0(P,  v)  =  R(P)  U(v)  , 
we  can  write  the  two  following  equations: 


V*R  -  -i-  gradpN  •  gradpR  -  IKN*  R  =  0 


where 

(Pv)  =N(P)  4  (v)  ,  f,  =  - gradp  fD. 

No 

and  Kis  an  arbitrary  constant  corresponding  to  the  separation  of  variables. 
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We  intend,  by  solving  these  equations,  to  analyze  the  influence  of  the 
sole  density  gradi<  nt  on  the  electromagnetic  properties  of  the  plasma.  As  far 
as  the  attack  to  the  system  (a.  10)  in  its  more  general  form  is  concerned,  or, 
at  least,  with  the  only  hypothesis  of  time-independence,  we  may  introduce; 

i.  some  possible  reductions  of  the  system  to  special  forms  which 
one  can  hope  to  integrate  in  closed  form, 

ii.  we  may  look  for  the  evaluation  of  eigenfunctions  and  eigenvalues 
of  the  opt  rator; 


a 

_ r 

/  MV- 

ov 

- 

\  J  rnv  <3  v  J 

•J 

to  integrate  the  system  by  means  of  expansion  processes  provided  that  the 
existence  of  a  closed  spectrum  for  J  can  be  proven. 

iii.  a  typical  iterative  process,  on  expanding  fy ,  fj,  formally,  in 
series  of  some  com  enient  parameters, 

iv.  iinally,  we  may  integrate  the  system  numerically,  at  least  in  the 
simple  cases  mentioned  above. 

b.  We  intend  also,  to  analyze  as  accurately  as  possible  the  validity  of  the 
approximations  made  in  the  spherical  harmonic  expansion  used  for  the  electronic 
distribution.  To  this  end,  we  will  eventually  try  to  consider  also  the  effect  of 

i 

the  second  order  anisotropy.  Furthermore,  within  the  limits  of  elastic  collisions, 
we  will  analyze  more  precisely,  the  particular  kind  of  electron-molecule  inter¬ 
action  which  affects  the  form  of  functions  vj  (v)  and  v^(v)  of  Section  l. 
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c.  As  regards  the  extension  of  the  theory,  vve  may  take  into  account 
more  than  one  kind  of  spherical  molecule  The  more  useful  extension,  however, 
concerns  the  possibility  of  taking  into  account  also  anelastic  collision  processes 
(excitation  and  ion. cation)  and,  at  ’east  partially,  and  in  the  limit  of  the 
Lorentnian  gas  model,  tee  remaii  ing  complex  t  f  ctron-moi  cule  interaction 
phenomena,  w  ii  n  are  pn  sent  i  1  the  plasma  under  various  physical  conditions 
(recombination,  molecular  attu -timc-nt,  et<  .). 
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LIST  OF  SYMBOLS  -  CHAPTER  II 


-  light  velocity  in  vacuum 
=  Boltzm.mi  1  s  constant 

=  elect  ron  mass 
=  c  le  1. 1  ron  c  ha  r ge 

-  distribution  functions 
=  eL'ctron  veloc  itv 

=  modulus  of  !  Irct  i'u:.  velocity  (v  =  |  j  ) 
=  molecular  mass 


molec  ul.ir  tempo rature 

modulus  of  La  ■  molccul.tr  \olocit\  (V  -  V|  ) 

molecular  velocity 

molecular  lle-a:.  \<  lot  it  y 

molecular  meal,  square  Srloiily 

elect ric  field 

in.tgm  tit.  field 


=  Boltzmann  differential  opt  rat  or  lor  electrons: 


Di  = 


Ui 

dt 


.  i'nd  ^  f 


qQ 


( i-  -  —  1  A  H  ) 


=  Boltzmann  ditle rential  operator  for  molecules: 


OF  ,  * 

DF  = - -  V  *  grad,-  i  r - 

ft  _  h  P  M 


gradv  F 


X  =  external  forces  on  molecules 


•  grad  f 
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CHAPTER  III 


MACROSCOPIC  PROPERTIES 

1.  INTRODUCTION 

For  this  chapter  the  preliminary  results  of  the  study  of  the  electro¬ 
magnetic  properties  of  a  nonuniform,  mac roscopic ally  neutral  plasma  are 
reported. 

The  line  of  approach  is  a  macroscopic  one,  based  on  the  thermodynamics 
of  irreversible  processes.  A  three-fluid  model  of  the  plasma  is  assumed  and, 
in  this  first  stage,  the  plasma  is  considered  to  be  at  rest  in  a  suitable  reference 
fiame,  the  components  are  in  mutual  thermal  equilibrium  and  imposed  magnetic 
fields  are  absent.  The  nonuniformities  are  those  connected  with  the  presence 
of  gradients  of  state  parameters  such  as  temperature,  pressure  and  concentra¬ 
tions. 

The  presence  of  these  gradients  induces  a  polarization  in  the  medium 
even  in  the  absence  of  external  electric  fields.  It  will  be  shown  that,  in  the 
stationary  state,  the  induced  polarization  depends  on  only  two  independent  gradients, 
which,  for  convenience,  are  taken  to  be  those  of  the  temperature  and  the  pressure. 
An  expression  for  the  induced  polarization  is  derived,  which  is  valid  for  a 
most  general  case  and  which  is  later  simplified  to  the  case  of  the  imperfect 
Lorentz  gas  (  i.  e.  weakly  ionized  plasma  for  which  the  mass  ratio  |  between  the 
negative  charge  carriers  and  neutral  molecules  is  much  smaller  than  one)  and  the 
perfect  Lorentz  gas  (  i.e.  weakly  ionized  plasma  with  (,  -  0)  . 
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The  phenomenological  coefficients  relating  the  induced  polarization 
to  the  temperature  and  pressure  gradients  are  first  expressed  in  terms  of 
thermodynamic  properties  of  the  plasma.  It  is  then  shown  how  these  coef¬ 
ficients  can  be  expressed  only  in  terms  of  simple  binary  diffusion  coefficients 
and  thermal  diffusion  coefficients  pertinent  to  the  plasma  constituents.  This 
last  step  furnishes  the  linking  element  between  the  present  macroscopic  approach 
and  the  more  refined,  but  usually  less  general,  statistical  approach.  Indeed 
it  indicates  how  the  subject  phenomenological  coefficients  can  be  rigorously 
evaluated,  in  tern's  of  the  "microscopic"  characteristics  of  the  plasma  consti¬ 
tuents,  once  the  statistical  approach  has  gone  far  enough  to  provide  workable 
results.  In  the  meantime,  however,  the  fact  that  one  needs  to  know  only  the 
binary  transport  coefficients  provides  a  possibility  of  readily  performing 
order  of  magnitude  analysis  by  utilizing  the  available  experimental  data  on 
electron-ions,  electron-molecules  and  ion -molecules  collisions. 

The  results  of  the  present  stage  of  the  analysis  already  open  up  two 
possible  ways  of  utilization,  both  of  noticeable  interest. 

In  the  first  place,  as  mentioned  before,  several  order -of-magnitude 
comparative  analyses  can  be  carried  out.  For  instance,  one  can  evaluate  the 
relative  importance  of  the  different  -inhomogeneities  with  respect  to  the 
polarizations  they  are  able  to  induce  (  e.  g.  ,  is  a  temperature  or  a  pressure 
gradient  more  important  to  this  effect?)  Also,  one  can  compare  the  polari¬ 
zation  induced  by  the  inhomogeneities  with  that  induced  by  an  external  electric 
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field  and  thus  determine  the  ranges  within  which  the  former  one  can  be 
neglected  and  the  medium  treated  as  homogeneous. 

The  second  application  of  equal  practical  interest  follows  from  the 
consideration  of  the  electromagnetic  field  which  is  created  by  the  induced 
polarization  (which,  in  principle,  can  be  computed  by  evaluating  the  Hertz 

jTTN 

potential  whos*  source  is  the  induced  polarization)  .  When  we  induce  into  the 
plasma  a  polarization  by  means  of  controlled  suitable  inhomogeneities,  part 
of  the  energy  furnished  to  maintain  the  gradients  is  spent  to  "create"  the 
induced  electromagnetic  field.  This  is  nothing  but  the  underlying  idea  of 
energy  conversion  systems  and  the  present  results  can  be  used  to  investigate 
both  the  feasibility  and  efficiency  of  a  number  of  such  systems. 

The  steps  of  the  analysis  herein  reported  are  as  follows: 

In  Section  2  the  irreversible  thermodynamic  description  of  the 
system  is  performed.  The  extensive  and  intensive  state  parameters  and 
the  pertinent  mass  and  energy  fluxes  for  a  mixture  of  three  lluids,  of  which 
two  have  negative  and  positive  charges,  are  defined.  The  'kinetic  relations 
between  fluxes  and  generalized  forces  are  established  and,  through  suitable 
use  of  the  basic  theorems  of  the  thermodynamics  of  irreversible  processes, 
the  general  expression  for  the  mass  fluxes  in  terms  of  pressure,  temperature, 
the  electric  potential  (6  )  gradients  and  thermodynamic  properties  of  the 
plasma  is  arrived  at. 
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In  Section  3  the  general  expression  for  the  polarization 
is  obtained,  for  a  macroscopically  neutral  plasma,  in  the  form: 

grad  6  =  o  j  grad  T  +  n  2  grad  p 

where  the  «  ^  are  phenomenological  coefficients.  Their  expression  in  terms 
of  the  thermodynamic  properties  of  the  plasma  is  also  derived  in  this  section, 
both  tor  the  general  case  and  for  the  simplified  case  of  the  Lorentz  gas. 

In  the  last  section,  Section  4,  the  problem  of  relating  the 
coefficients  a  to  binary  transport  coeficient  is  considered.  By  suitable 
transformation  of  fluxes  and  affinities  it  is  shown  how  they  can  be  expressed 

in  terms  of  three  binary  diffusion  coefficients  D2  >  and  two  thermal 

1'  1' 

diffusion  coefficients,  D  ,  D,  which  refer  to  the  plasma  constituents. 
General  expressions  are  presented  but  no  attempt  is  yet  made  to  either 
simplify  them  or  evaluate  their  comparative  order  of  magnitude. 

The  concluding  section  briefly  describes  the  lines  along  which 
the  research  will  be  furthered 
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IRREVERSIBLE  THERMODYNAMICS  DESCRIPTION  OF  THE  MEDIUM 


2. 


Considering  a  plasma  formed  by  electrons  (  subscript  1  )  ,  ions  (  l)  and 

neutral  molecules  (  3)  in  neutral  thermal  equilibrium  (  i.  e.  T,  =  T2  =  Tj  =  T) 

one  can  write  the  mass  fluxes  for  the  three  species  and  the  energy  flux  as  a 

linear  combination  of  generalized  forces  (  affinities)  through  phenomenological 

coefficients.  The  plasma  is  assumed  to  be  at  rest  in  a  suitable  reference  frame) 

As  a  consequence  of  the  mass  conservation,  only  two  mass  fluxes  are 

independent  so  that,  by  accounting  for  the  Onsager's  principle  of  the  symmetry 

of  the  phenomenological  coefficients  (  L-  =  Lj-)  ;  one  can  write: 

2  '  7  J  v 

Jo  *  Y  °,Ji  * 


is  1 


oo 


nts  (  L-.  =  Lji)  ;  one 

o.y. 


(2.  1) 


2 

Ji  =  Y.  Lik  r>:k  -  X,  +  Qk  X.) 

k:  1  '  ' 


(  i  =  1.2) 


In  these  relations  Jq  is  the  flux  of  energy,  J-'s  are  the  mass  fluxes,  and  the  s 


are  the  generalized  forces: 
Xq  =  -  ^  grad  T 


X,  =  -e!  grad  «b  -  T  grad 


X2  -  -  e2  grad  <b  -  T  grad 


(H 

w 


(2.2) 


X3  =  -T  grad 
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with  p.  =  specific  chemical  potential  for  the  itd  species,  e  the  electric 

charge  per  unit  mass  and  <J>  the  electric  potential. 

The  quantities  Q.  (  i  =  1,2)  are  the  heats  of  transfer,  defined 

in  terms  of  the  coefficients  L;  .  as: 

J 

„  _  l-lO^Z  '  ^12^20 

Ul  ~  - T - 7~1 — 

J-l  1  t- 22  -  L,  2 

(  2.  3) 


Qz 


L1  1  L2Q  -  Lj  o  1-21 

L,  ,  L22  -  L,  2 


It  is  convenient  to  eliminate  from  Equations  (  2.  1)  the  explicit  appearance  of 

the  gradients  of  the  chemical  potentials  pk  in  favor  of  the  pressure  gradient. 

This  can  be  done  by  considering  the  p^  =  pk(T,  p,  Cj)  ,  where  Cj  is  the  mass 

th 

concentration  of  the  i  species,  and  by  taking  into  due  account  the  definition 
of  Gibb's  potential  (Reference  2)  and  the  Maxwell's  relations  (Reference  2)  . 


The  result  is: 
2  r 


J.  =  7  Lq  < 

1  kTi  *k 


114-  1 


( v3  -  vk)  grad  p  - 


j=  1 


0  ( Kk"  Kj) 

dCj 


1 


grad  c 


P’  T  ’  ck 


J 

(^.4) 


-  (  Qk-  \  +  hj)  t  grad  T  f  ‘  ei  grad 


J 


i  =  (1,2) 
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where  v-  and  h'  are  the  specific  volume  and  enthalpy  of  the  i1*1  component  and 
where  the  subscripts  indicate  the  variables  to  be  kept  constant  in  the  differen¬ 
tiation  process. 

Equations  (  2.4)  are  the  equations  needed  to  determine  the  polarization 
induced  by  the  i nhomogeneities  in  the  plasma,  as  detailed  in  the  next  section. 

3.  GENERAL  EXPRESSIONS  FOR  THE  POLARIZATION 

We  are  concerned  with  the  determination  of  the  polarization  (grad  6) 
induced  in  the  plasma  by  imposed  gradients.  We  deal,  therefore,  with  a 
stationary  state  (Reference  1)  wherein  all  the  mass  fluxes  will  be  zero  while 
the  energy  flux  J0  is  different  from  zero. 

Equations  (2.  4)  equated  to  zero,  plus  the  statements  of  mass  conserva¬ 
tion  and  macroscopic  electrical  neutrality  provide  a  set  of  four  equations  for 
the  six  quantities  grad  c.  (i  =  1,2,  3)  ,  grad  T,  grad  p  and  grad  d>  .  It  follows 

that  any  of  these  quantities  can  be  expressed  as  functions  of  two  others. 

By  choosing  the  two  independent  gradients  to  be  the  pressure  and  tempera¬ 
ture  gradients,  one  thus  obtains  the  following  general  expression  for  the  polari¬ 
zation: 


grad  6  =  a  j  grad  T  +  a  2  grad  p.  (3.1) 

The  coefficients  q,  and  a2  are  herein  referred  to  as  coefficients  of  polarizability. 
They  are  a  measure  of  the  polarization  induced  by  the  presence  of  gradients  of 
macroscopic  properties  of  the  plasma  and  are  expressed  in  terms  of  their 
thermodynamic  properties  by: 
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a  i 


with: 


(3.2) 


ai  =  9  (  Hq  '  1^3 )  /  9C, 
bi  =  a  (Ki  -  n3)  /  dc2 

™i  =  vi  -  v3  (  i  =  1.2) 

1 

"p  (  Oj  -  +  hj  ) 


(  3.3) 

* 


It  appears  that  the  a-,  b^,  are  completely  determined  once  the  pertinent 
"state  equations"  (Reference  2)  of  the  constituents  are  given.  The  quantities 
q^  depend,  in  addition,  upon  the  phenomenological  coefficients  through  the 
heats  of  transfer  Q.. 

Thus,  to  make  any  use  of  Equation  (  3.  2)  one  must  make  some  assumptions 
as  to  the  state  equations  of  the  plasma  constituents  and  must  determine,  either 


experimentally  or  by  means  of  statistical  mechanics,  the  phenomenological 
coefficients  L-  . 

M 
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The  expressions  for  a*.  b.  and  m-  are  herein  derived  on  the  assumption 

that  the  plasma  constituents  are  perfect  gases.  The  quantities  q.  will  be  dealt 
with  in  the  next  section. 

If  one  makes  the  perfect  gas  assumption,  the  following  relations 
hold  (Reference  2)  : 


v.  = 
1 


h.  = 
i 


m 


h- 


10 


RT 


Pci 


(  3.-1) 


wherein  (  T)  is  an  arbitrary  function  of  the  temperature,  m 
t  h 

mass  of  the  i  component,  m  is  the  molar  mass  of  the  mixture, 
specific  enthalpy,  and  &  the  molar  gas  constant. 

V/ith  these  expressions  one  has: 


is  the  molar 
hQ  a  reference 


to 
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In  the  case  ol'  an  imperfect  Lorentz  gas  (  i.e.  —I-  <  <  1;  N,  /N,  <  <  1  • 

in,  1  3  ’ 


N,/N3  <  <  1)  these  formulae  are  greatly  simplified  and  one  obtains: 


a,  NRT  m 

— ■*“  -  — 

e,  -  «  ,  m 


l 


b,  =  az  ~  0 

b3  ^  NRT  m 
e2  ~  £  i  m. 


(  i.  6) 


mi  = 


RT  mi 


pm  ct 


where  c  ,  and  e  ,  are  the  molar  electric  charges. 

The  expressions  for  the  coefficients  U]  and  n  >  are  found  by  substituting 

Equations  (  5.5  )  into  Equations  (  5.  . 

Another  quantity  of  practical  interest  is  the  coefficient  n3  -  u^/ci,  ,  giving 
i 

a  measure  of  the  relative  importance  of  pressure  and  temperature  gradients  with 
respect  to  the  induction  of  polarization. 

The  general  expression  for  q3  is: 


whe  r  e 


m,  +  m,  0 


q,  +  q3e 


±1_ 

e2 


(  i.  7) 


(  i.8) 


dL. 

e?  e. 


1 


1 

I 

I 

I 

I 


In  the  case  of  imperfect  Lorentz's  gas,  Equation  (  5.7)  reduces  to: 


RT 

m. 

€  ^ 

mz 

i 

„  +.  <  z  Inz 

u3  = 

pm 

[  c. 

«  1 

ml  _ 

i 

Lq'  ^  77  TTTj 

-  Nrn 


c  2ni2 


Cj  m 


_\ '  /  U 

')  V 


1 


(  5.  10) 


The  formulae  so  far  developed  are  all  that  one  can  do  without  specifying  the 
nature  of  the  phenomenological  coefficients  L3|.  Any  order  of  magnitude  analysis 
can  be  furthered  only  after  having  found  suitable  expressions  for  these  coefficients, 
or,  what  amounts  to  the  same,  for  the  quantities  q^. 

This  is  considered  in  the  next  section,  where  it  will  be  shown  how  the 
L •  's  can  all  be  expressed  in  terms  of  binary  transport  coefficients. 

j 

4.  EXPRESSIONS  OF  THE  PERTINENT  PHENOMENOLOGICAL  COEFFICIENTS 
IN  TERMS  OF  BINARY  TRANSPORT  COEFFICIENT 

To  express  the  coefficients  in  terms  of  binary  transport  coefficients  the 
following  consideration  is  essential. 

In  the  subject  case  we  have  two  independent  mass  tl  lxes  and  one  energy 

flux.  The  number  of  independent  phenomenological  coefficients  is  therefore, 

accounting  for  the  three  Onsager's  relations  between  the  six  c ros s -coetf ic ionts , 

six.  It  follows  that  the  behavior  of  the  system  is  completely  characterized  by 

only  six  independent  coefficients.  A  set  of  six  independent  coefficients  is  that 

of  the  thjree  binary  diffusion  coefficients  D,  ,  ,  D13,  D,3  (  i.e.  the  coefficients  of 

diffusion  for  electrons  -ions,  electrons -molecules  ,  and  ions-molecules  mixtures) 

T  T 

plus  the  two  binar  y  t  he  r  mal  diffusion  coefficients  D]  ,  D2  (i.e.  the  coefficients 
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of  thermal  transport  due  to  diffusion  of  the  electron-gas  and  the  ion-gas  into 
the  neutral  gas,  respectively)  and  the  heat  conduction  coefficient. 

It  is  then  natural  to  think  that  a  suitable  linear  transformation  of  fluxes 
and  forces  will  make  it  possible  to  express  everything  in  terms  of  the  above 
binary  transport  coefficients. 

This  is  indeed  so.  In  fact,  by  defining  the  new  fluxes  Jk  and  forces  X^. 
through  the  relations: 

j/  =  J,  (i  =  1,2.3) 

Jo  =  Jo  ‘  hl  J>  ‘  hlJz  '  Jj  (4'  1} 

x;  =  x.  +  hiXo 
x’  =  x 

O  O 

We  find  (  see  Reference  2)  the  following  relations  between  the  coefficients 

t 

]_■  and  the  new  coefficients  L-.  . 


Hi  - 

Lil 

(i.j  =  1.2,  1) 

Loo  - 

Loo 

(4.  2) 

L>o  ^ 

-Ln  (h, 

-  h3)  -  L,  j,  (  h2-h3) 

*  ^1  o 

1 

f-’ZO  = 

-L2.  (h, 

-hs)  -  -  h3) 

4'  I-1’  o 

1 

These  new  coefficients 

are  just  those  used 

in  Reference 

3  and  it  is 

therein  shown  how  they  are  related  to  the  binary  diffusion  coefficients  Di- .  ,  to 
the  thermal  diffusion  coefficients  D-  (  i  =  1,2)  and  to  the  thermal  conductivity 


coefficient  h'. 


By  performing  the  necessary  substitutions,  the  ultimate  goal  of 
expressing  the  q-'s  in  terms  of  binary  transport  coefficients  is  thus  achieved 


and  one  obtains: 


q!  = 


-  AT  *  n^ni  n2™i  m2  Jnjm3  d;j  D23  +  (p=n2m2)  D,  ,  Da  f 

.  q  r  T  T"1  T  p  2 

-  m,  (  p  -  n,  m,  )  Dj  2  D.  3  D,  +  D,  +  (  Dl  )  n‘n2n3  m2m3  n,  m,  D21  D. 


-m3  (  p  -  n3  m}  )  Da  DJ;  -  m2  (  p  -  n,  m, )  Da  D21 


and 


= 


AT 


n  n,  n,m,  m> 


[  n3  mi  D1  3  D23 


-  ni2  (  p  -  n2m2)  Dj  2  Da  + 


-m,  (  p  -  n,  irij 


D14n,jJ  ID 


T  T 

l  +  D2 


(  D, 


•J 


•  1  2 

n,m  D;2D32  -  m3(p-n3m3) 

Di  3  Dj2 

•m,  (  p*n:  m,  ) 

D.jDjJj  > 

where  A  is  given  by: 

p* 

✓ 

A  = 

P"P  <  n,  D,j  +  n2  D,  3  +  n}  Dj  2) 

~  1 

•  < 

2 

1  n  nj  n}  m :  m3 

n2m2D12DJ2  + 

‘1  n2ml  m2 


•  nij  (  p  -  iij  m3 )  D;  j  Dj2  -  nij  (p-njtrij)  Dj  3  D,  2^J  n’lij 
> 

n3  ,n3  Di  3  D>3  -  m2(  p-  n2  ni2 )  D,  2  D23  -  m,  (  p-  n3  m,  )  D,  .  Dj  3 

J^n2  m2  Dj  >  Dj2 -m3  (  p-n3  m3)  D  3  D3 ,  -m,  (  p -n,  m,  ) 


n  rij  n3  nij  m3 


•  Dj  3  Dj  2  j  n2n2n3m2m3  n,  mf  D2I  D} ,  -m,  (p  -n3m3)  DaDj,  + 


-m2  (  p  -  n,  m,)  D21  +  n2  nj  n2  m,  m. 


n3ni3Di3D23  -  m2  (  P*  n2m2)  • 
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Dj  2  Dy  -  m,  (  p  -  n,  m,  )  D;  3  D[  2  n*  n2  n3  m2  n 


ni  m,  D21  D3i  + 


-m3  (  p  -  n3  m3  )  Da  D31  -  m2(  p  -  n2  m2)  Dv,  D 
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At  this  point  the  new  polar ization  coefficients  result  expressed  as  functions 
T 

of  D-  and  D.  only.  One  could  carry  the  analysis  further  and  relate  them  to  the 

\>  1 

appropriate  integrals  of  the  distribution  functions  for  the  three  constituents  thus 
establishing  the  link  between  the  present  macroscopic  approach  and  the  statistical 
one.  We  don't  pursue  this  aspect  any  further  since  all  the  required  passages 
and  formulae  are  developed  extensively  in  Reference  1,  (part  L,  Chapter  7)  to 
which  the  reader  is  referred  for  details. 


5.  CONCLUDING  REMARKS 

The  general  expression  for  the  polarization  induced  by  inhomogeneities 
has  been  herein  derived.  Two  new  "polarizability"  coefficients  nj  and  n2  have 
been  defined  which  are  a  measure  of  these  effects  and  which  have  been  expressed 
in  terms  of  the  thermodynamic  stale  equations  of  the  plasma  constituents  and  of 
binary  transport  coefficients. 

In  the  subsequent  stage  of  the  research,  an  order  of  magnitude  analysis  will 
be  performed  to  determine  both  the  absolute  and  relative  importance  of  temperature 
and  pressure  gradients  with  respect  to  the  induction  of  polarization  in  a  plasma. 

This  can  be  done  by  a  judicious  choice  of  the  values  of  the  binary  transport 
coefficients,  as  they  are  presently  available  in  the  literature  on  the  assumption 
that  the  constituent  gasses  obey  perfect-gas  state  equations. 
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The  present  method  of  approach,  furthermore,  will  be  extended  to 


include  the  following  effects: 

i)  presence  of  an  imposed  magnetic  field  (i.e.  study  of  a  plasma 
medium  which  is  both  nonhomogeneous  and  anisotropic) 

ii)  plasma  constituents  not  thermally  in  equilibrium  (  i.  e.  electron,  ion 
and  neutral  gasses  with  different  temperatures) 

iii)  nonuniform  convective  motion  of  the  plasma. 
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SYMBOLS  -  CHAPTER  III 


mass  concentration  of  i  component 
T 

D  coefficient  of  thermal  diffusion 

i 

D..  coefficient  of  binary  molecular  diffusion  between  the  species  i  and  j 
electric  charge  per  unit  mass  of  the  i  component 
E  external  electric  field 

g  Gibbs'  specific  potential 

h  specific  enthalpy 

J-  mass  flow  of  i  component  referred  to  the  baricentcr  of  the  mixture 

Jq  ene  rgy  flow 

€  rm  N. 

m  molar  mass  of  the  mixture  =  - - 

N 

N 

n  number  density  of  the  mixture  = 

n.  number  density  of  i  component 

N.  molar  concentration  of  i  component 

1 

N  molar  concentration  of  the  mixture 

p  pressure 

P  polarization 

s  specific  entropy 

T  temperature 

u  specific  internal  energy 

V  specific  volume 

dp  a2  new  polarizability  coefficients 
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SYMBOLS  (  Continued) 


f  molar  electric  charge 

thermal  conductivity  coefficient 
p.  electrochemical  potential 
p  density 
<t>  electric  potential 

u)  angular  wave  frequency 


CHAPTER  IV 


PROPAGATION  OF  ELECTROMAGNETIC  WAVES  IN  THE  SHOCK  TUBE 
1 .  INTRODPCTION 

The  present  part  of  the  theoretical  work  on  the  propagation  ot  an 
elect romagnetic  field  in  a  plasma  is  connected  with  the  experimental  program 
of  measurement  of  the  electrical  properties  of  the  plasma  in  the  shock  tube. 

The  test  section  of  the  shock  tube  is  used  as  a  section  of  wave  guide  in 
which  the  phase  velocity  and  the  attenuation  of  the  selected  mode  of  propagation 
will  be  measured.  The  propagation  constant  of  the  electromagnetic  wave  depends 
upon  both  the  axial  and  radial  distribution  of  the  electric  properties  of  the  plasma. 

The  axial  gradient  of  electron  density  is  controlled  by  the  temperature  distribution 
in  the  driven  section  of  the  shock  tube.  On  the  other  hand  we  may  expect  a  radial 
nonuniform  distribution  in  the  electron  density  due  to  the  growing  thickness  of  the 
boundary  layer  behind  the  shock  wave.  The  rate  of  increase  of  the  boundary  layer 
thickness  corresponds  to  a  small  change  of  the  radial  distribution  inside  a  wavelength 
of  the  field  distr ibution  in  the  wave  guide.  Consequently  the  effects  of  radial  and 
axial  nonunil'orm  distributions  can  be  analyzed  independently.  .  . 

.  The  local  value  of  the  phase  velocity  along  the  wave  guide  depends  upon 
some  average  value  of  radial  electron  distribution.  This  dependence  must  be 
analyzed  in  order  to  determine  the  error  which  may  be  introduced  in  the  correlation 
of  the  experimental  results  with  a  propagation  model  in  which  the  radial  distribution 
is  not  taken  into  account. 


?  5 


r—t 


In  order  to  discuss  this  effect  a  theoretical  analysis  of  the  propagation 
properties  is  being  conducted  assuming  a  circular  wave  guide  of  radius  r^  in 
which  the  plasma  is  contained  in  a  cylinder  of  radius  r,  <  r_>.  The  values  of 
r;  -  r,  are  selected  according  to  the  expected  thickness  of  the  layer  in  which 
the  electron  density  decays  rapidly  in  the  nc ighbo rhood  of  the  wall.  The  step 
function  distribution  is  selected  in  order  to  obtain  the  maximum  possible  effect 
of  the  boundary  layer  on  the  axial  propagation  of  the  electromagnetic  wave.  This 
problem  is  analyzed  in  Section  i  of  this  chapter.  In  the  final  stage  of  the  ex¬ 
perimental  program  wo  arc  planning  to  measure  the  propagation  properties  in 
the  wave  guide  also  with  an  axial  constant  and  uniform  magnetic  field  in  the 
test  section  of  the  shock  tube.  These  measurements  will  provide  an  additional 
set  of  data  on  the  electromagnetic  properties  of  the  plasma.  For  this  reason 
a  theoretical  analysis  is  being  conducted  in  order  to  calculate  the  propagation 
properties  and  in  particular  to  discuss  the  effect  of  the  magnetic  field  on  the 
lowest  modes  of  propagation  of  the  waveguides  containing  the  ionized  gas.  The 
preliminary  results  of  these  calculations  are  contained  in  the  last  section  of  this 
c  hapte  r . 

i.  BASIC  EQUATIONS 

The  mathematical  model  of  Maxwell  equations  together  with  an  electron 
momentum  equation  is  developed  first  in  a  general  way  and  then  specialized  to 
the  specific  problem  treated.  The  plasma  is  assumed  to  be  defined  by  a  constant 
plasma  frequency  and  a  constant  collision  frequency.  The  system  of  equations 
upon  which  the  analysis  is  based  is  as  follows: 


I 


^6 


V  x  E 

=  -  Hp 

(  1) 

V  X  IT 

=  «  E  -  e  n  w 

(  £) 

• 

m  w  ~ 

-  v  m  w  -  e  E  -  u.  e  w  X  H 
^  o 

(  3) 

where  E,  H,  are  respectively  electric  and  magnetic  oscillatory  fields  and  H0 
is  the  axial  applied  (  steady)  magnetic  field,  w  is  electron  velocity,  v  is 
collision  frequency,  m  is  electron  mass,  (  -e)  is  electron  charge,  and  e  ,  p 
are  dielectric  permeability  and  magnetic  permeability,  respectively,  for  a 
vacuum.  The  electrons  are  assumed  to  have  significant  collisions  with  neutral 
particles  only,  and  v  is  the  frequency  for  this  type  of  collision. 

Under  harmonic  time-dependence  e*1^  *•  the  equations  become: 
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The  curl  operator  is  partitioned  in  the  following  manner: 


i 

j 

k 

i 

j 

k 

Vx  = 

d_ 

dx 

a 

ay 

0 

+ 

0 

0 

a 

3z 

(  ) 

(  ) 

(  ) 

(  ) 

(  ) 

(  ) 

Vs  X  (  )  +  Vz  X  (  ) 

where  the  partitioning  is  illustrated  in  Cartesian  coordinates  for  simplicity, 
z  still  being  considered  as  axial. 

Applying  this  partitioning  to  the  system  (1  )  ,  ( 5  )  ,(  6)  ,  the  following 
system  results: 


VsXEs  *-  i«nHz 

V  g  X  Ez  ,  Vz  X  Es  =  iu  |xlls 

Vs  X  Hg  =  -iu.  «  E  -  enwz 

vs  X  Hz  +  V  X  IIS  =  -iu  <  Es-enws 
e 

(  v  -  i  u  )  wz  =  -  -  Ez 

(v  -  iu.  )  ws  =  -  ^  Es  x  u  Hiz  X  wg 


(9) 

(  10) 

(  ID 

(  12) 

(  13) 

(  H) 


In  the  above  system,  Equations  (  9)  ,  (11)  ,  (13)  are  scalar  single 
equations  and  (  10)  ,  (  12)  ,  (  14)  are  vector  equation  pairs  in  the  (  x,  y)  or 
(  r,  @)  plane.  The  last  equation  (  14)  can  be  solved  for  the  two  components, 
wr  and  Wg,  of  wg  : 
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I 


I 

I 

i 

] 

i 


{  1  5) 


or 

where 


(lb) 


T  I, 

\  *  u  H 


F rotn  (14) 


(  17) 


w 


z 


-e 

m  (  v  -  i  at  ) 


E„ 


(  18) 


Substituting  from  (  1 6)  and  (  18)  into  (  11)  ,  (12)  for  w^,  \v  ,  there  results 


V 

s 


X  1 1, 


i  w  «  E  + 
z 


e2n  E 

IT}  (  V  -  io  ) 


i  w 


E, 


+ 


E 


z 


R  E 


z 


who  re 


—  is(plasma  frequency) 
r  m 


V  X  H  +  V  X  H 

S  Z  /.  s 


i  u  €  Eg 


- 1 

T 


E 

s 


=  R  E 


s 


where  the  scalar  R  and  the  matrix  R  are  given  by 


R 


v  ■  iu 


R 


-  i  u>  €  I  + 


-l 

T 


with  I  representing  the  2X2  unit  matrix. 


(  19) 

(  ZO  ) 

(  Zl) 

(  Z2) 


39 


Now  solving  formally,  (  10)  for  II,  and  (  20)  for  Eg,  we  find 


H 


1 

i  u  p 


(  V  X  E  +  V  X  E  ) 


(  23) 


Eg  =  R  (  Vs  X  Hz 


V,  X  "s  > 


(  24) 


The  objective  here  is  to  express  the  transverse  fields,  Hg  ,  Eg  entirely  in  terms 
of  the  axial  components  E  ,  H  .  This  is  accomplished  by  substituting  (  23)  into 

Z  Z 

(  24)  and  (  24)  into  (  23)  . 


H  = 


s  1U  fl 


(  V  X  Ez  +  V  ?  X  (ll‘‘  X  Hz  +  R'1  Vz  X  Hg  )  )  (25) 


E  =  R*1  (  V  X  H  +  V  X  (  V  X  E  +  V  X  E  )  \ 

s  \s  ziupz  s  z  i.  s  / 


(26) 


But  it  can  be  easily  shown  that 


V  X  Vz  X  Es 


d'E, 


Oz 


(  27) 


and 


V  X  (  R'1  V  X  H  ) 
z  z  s 


-i  d  H 

R  _ !_ 


dz. “ 


(  28) 


In  the  cavity  mode  problem,  only  the  first  axial  harmonic  is  desired.  Thus 
for  the  cavity  problem  the  z-dependence  of  Ezand  I^is  given  by 


I 


Ez(  r,  6) 
Hz(  r,  6  ) 


TT  Z 
Sill  — 

z 

o 


1TZ 

COS  - 

z 


o 


z 


(  29) 

(  30) 


so  that  9  (  )  ~  ^  1 _ 1 


where  z  is  the  axial  length  of  the  cavity.  In  the  wavequide  problem  the 
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z-dependence  of  Ezand  I ^ i s  given  by 


E5,  (  r ,  6)  e 


IV  z 


Hz(  r,  0)  e 


iVz 


{  32) 
(  33) 


so  that 


li 

f)z* 


(  34) 


Hence,  in  what  follows  d’/oz”  may  be  replaced  by  either  -  it2/zq  or  -'y1 
and  no  confusion  should  result.  This,  together  with  (  2.1 )  and  (28) 
permits  solution  of  (  25)  and  (  26)  for  E  ,  H  in  terms  of  E  ,  H  . 

^  S  2  ^ 


E 


t“  R  V 

I  -  ~T— 

,  V  w  11  J 


r'1  (  V  X  H  +  -1 —  V  XV  X  E 

‘s  z  lupz  s  z 


(  35) 


i  it“R  * 1 


- 1 


“.-.'•‘Tt 


1 


»  \ 


zolu  V-. 


iu  p 


Vc  X  E,  +  V.  X  (  R_1  V  X  H 
“  z  z  s  z 


(36) 


This  expresses  Eg  and  Hg  in  the  desired  form. 

For  convenience  in  subsequent  expressions,  the  following  notation  is 
introduced 


R'1 


(  37) 


(  38) 


Then,  if  E  and  II  from  (  35)  and  (36)  are  substituted  into  (9)  and  (  19)  , 

’  s  s 

using  the  n  ,  (3  ,  V  >  b  notation  of  (  37)  and  (  3  8)  ,  the  following  two  coupled 
equations  for  H  (henceforth  called  simply  H)  and  Ez  (called  E)  are  obtained: 
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Letting  the  0-dependence  be  given  by  e1"  *  in  all  cases,  and  using  the 
operator  notation 


—  ±  (  r  r) 
ni  r  dr  dr 


m 


(41) 


the  equations  (  39)  and  (40)  reduce  to 


Rru  (J.E  +  6L  E  1  —  (  \  n  t-  6  |3  )  L,  ,  H  -  0 

ni  z  1  '  m 


4  -) 


iu»  p  H  +  (  60  -i,  p  )  L  II  -  —  (  ^  -  6(3  )  L  E  -  0  (43) 

m  i  w  p  11 


Eliminating  H  by  cross  operating  on  these  equations,  we  find 


a  L  E  +  b  L  E  +  cF.  =  0 

m  m 


(44) 


with 
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This  has  the  solutions  Jm(  rk,  )  ,  J  _(  rk^)  ,  Ym(  rk,  )  ,  Ym(  rk,) 
where 

k. 
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(  17) 


(48) 


Incidentally,  H  can  be  shown  to  satisfy  the  same  equation  (14)  and  hence  to 
have  solutions  composed  of  the  same  functions. 

The  combinations  of  a  ,  |i ,  \  ,  6  required  are  as  follows  for  the  case 
of  zero  collision  frequency  (  v  ~  0  )  . 
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Note  that  p  and  (  y  a  4 


6  p  )  vanish  with  the  magnetic  field  (  i.  c.  as  u>  ^ 


->  0)  . 
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3 .  PROPAGATION  OF  AN  ELECTROMAGNETIC  WAVE  IN  A  CIRCULA R 
WAVEGUIDE  WITH  A  PIECEWISE  CONSTANT  RADIAL  ELECTRON 
DENSITY  DISTRIBUTION 


The  sketch  shows  the  radial  va  ria  tion  of  electron  density  considered  in  the 
wave  propagation  analysis.  As  noted  earlier,  this  is  intended  to  simulate  the 
rapid  radial  decay  of  electron  density  in  the  boundary  layer.  For  this  problem, 
no  externally  applied  steady  magnetic  field  is  considered.  The  uncoupling  realized 
in  the  cavity  mode  problem  by  the  disappearance  of  applied  magnetic  field  does 
not  carry  through  to  the  present  problem  (except  for  0-indc  pendent  waves)  due 
to  the  radial  nonuniformity  of  the  plasma  electron  density. 

However,  with  the  disappearance  of  u>H  ,  the  matrices  T,  R,  etc.,  degenerate 

to  simple  scaler  quantities.  Thus 
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while  a 
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In  the  above  expressions  (50)  through  (53)  the  value  of  wp2  is  to  be  taken 
<is  w^2  in  the  plasma  and  zero  in  the  vacuum  annulus. 

With  the  vanishing  of  and  y,  equations  (-42)  and  (43)uncouple,  and  both 
E  and  H  satisfy  a  reduced  wave  equation  of  the  form 
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u  representing  E  or  H. 

These  have  the  solutions 


u  = 


Jm  (l0r  .  if) 

\  u 2  +  *  ko  / 
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(55) 


In  the  outer  vacuum  annulus,  both  solutions  (55)  must  be  used  (for  both 
E  and  H)  with  idp  =  0  .  In  the  inner  plasma-filled  cylinder,  only  the  Jm  contribu¬ 
tions  apply,  from  boundedness  considerations.  This  leaves  three  undetermined 

V 

coefficients  for  each  of  E  and  H.  Correspondingly,  there  are  two  boundary 
conditions: 


E(  r2  )  =  0 

E(0)(ri)=° 


(56) 
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(57) 


and  four  interface  or  jump  conditions  at  r  =  r,  . 

[e]  =  0 
0 

w 

In  (5  ()  the  symbol  ^  j  means  jump  in  a  quantity  across  r  =  r,  . 

e  jump  conditions  (57)  on  E^  and  require  expression  of  these 

quantities  in  terms  of  E  and  H.  From  equations  (35)  and  (36)  and  the  scalar 
nature  of  R  (50),  the  expressions  for  E(0)  and  H(Q)  can  be  written  as  follows; 


=  -  6q  i1!  +  JLi  1  92E 


'(») 


3  r 
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H  =  -  6  +  6  n  _l  d*H 

iup  rir  iwp  r  az80 


(59) 


Now  equations  (58)  and  (59)  show  that,  in  absence  of  9 -dependence,  the 
problem  uncouples  into  TE  and  TM  modes,  while  the  jump-condition  require¬ 
ments  (57)  prevent  this  uncoupling  in  the  presence  of  9 -dependence. 

The  jump  and  boundary  conditions  now  lead  to  the  dctcrminantal  equation 
which  relates  propagation  velocity  7,  to  wave  frequency  w.  and  the  plasma 
frequency  wp  and  the  configuration  r,  /r;  .  After  some  algebraic  manipulation 
this  dete rrninuntal  equation  can  be  put  in  the  following  form;  r° 

i  /  / 
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l  \  m 


(60) 
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where  the  parameters  used  are  as  follows: 
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The  equation  (60)  constitutes  a  transcendental  relation  between  u,u>p,v,V» 
and  the  configuration  ^r, /r2 ) »  This  equation  is  currently  being  solved  numerically 
on  the  GASL  Bendix  G-15  computer. 
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4.  RESONANT  MODES  IN  PLASMA-FILLED  CAVITY 

In  this  section,  the  effect  of  an  external  axial  magnetic  field  is  analyzed 
assuming  a  standing  wave  distribution  of  the  electromagnetic  wave  inside  the 
wave  guide.  In  this  case,  we  may  study  the  resonant  modes  of  a  plasma-filled 
cavity  of  length  zQ  and  radius  rQ,  under  the  ext.  rnally  applied  magnetic  field 
H  The  minimum  value  of  the  length  zQ  of  the  cavity"  corresponds  to  half  a 


wavelength  of  the  standing  wave  inside  the  wave  guide. 

The  Y  solutions  must  be  ruled  out  from  considerations  of  boundedness 

at  the  origin.  Thus,  for  this  problem,  E  is  given  by  , 

i(-  wt  +  m#  ) 


IT  Z 


E  =  (A  Jm  (  r  k,  )  t  B  Jrn  (r  k,  )  J  sin_£-  e 


(6  2) 


The  remaining  two  boundary  conditions,  from  which  A  and  B  (the  mode  shape) 
and  the  resonant  frequencies  are  obtained  are: 

E  (  r0  )  =  0 

and  E^(  r°  ^  =  ® 

The  first,  (63),  follows  simply  from  (63): 

A  J™  (r0k,)  B  Jm  (r0k,  )  =  0 
The  second,  (64),  requires  the  use  of  (36)  for  Es,  which  in  turn  requires  a 
knowledge  of  H.  Using  (43)  and  (43),  H  is  expressible  in  terms  of  E  and 
operators  thereupon. 

When  this  process  is  carried  out,  and  when  the  vanishing  of  E  itself  on 
the  wall  r  =  rQ  is  accounted  for,  the  second  homogeneous  equation  on  A  and  B 

becomes: 


(63) 

(64) 

(65) 


t 
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L  H  1  w  K 

Equations  (65)  and  (66)  provide  the  determinantal  relation  for  a  nontrivial 


solution  vector  (A,  B)  from  which  the  modal  frequencies  are  computed.  It  can 
easily  be  shown  that,  for  =0  (zero  applied  magnetic  field)  k,  and  k2  become 
equal,  invalidating  the  solution  form  (  6^) .  However,  examination  of  (4Z)  and 
(43),  (noting  that  and  y  vanish  with  u>H)  shows  that  H  and  E  decouple  in  this 
case,  permitting  decomposition  into  the  conventional  TM  (H=  0)  and  TE(E=  0) 
modes.  The  natural  modes  and  frequencies  for  this  case  (when  v  =  0)  are  classical 
provided  w2  in  the  classical  problem  is  replaced  by  (w2  -  wp2  ).  That  is 


w(n)  “  w(n)2classical  +  wp2  (67) 

t  h 

for  the  n  mode. 

The  determinantal  relation,  corresponding  to  equations  (65)  and  (66)  is 

currently  being  solved  numerically  in  the  G ASL  Bendix  G-15  Computer.  At  a 

set  of  values  of  wH>  wp  and  geometry  rQ/zQ,  trial  values  of  frequency  w  are 

inserted  until  a  zero-crossing  is  obtained  for  the  value  of  the  determinant 
( 

Great  care  must  be  taken  in  the  interpretation  of  the  solutions  so  obtained,  since 
some  formal  solutions  exist  (vanishing  of  the  determinant)  which  correspond  to 
trivial  fields.  Such  is  the  case,  for  example,  when  either  kj  or  k,  vanishes; 
then  it  is  possible  to  have  a  trivial  field  even  though  the  vector  (A,  B)  is  nontrivial. 
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CHAPTER  V 


EXPERIMENTAL  FACILITY 

L  INTRODUCTION 

In  this  chapter  we  describe  the  facility  which  will  be  used  for  the 
measurement  of  the  electromagnetic  properties  of  a  nonuniform  plasma. 

The  first  part  contains  a  description  of  the  shock  tube  and  the  results  of 
calculations  concerning  the  state  of  the  test  plasma.  The  second  part  con- 

tains  a  description  of  the  test  section  and  the  discussion  of  the  electromagnetic 
measurements. 


2.  THE  SHOCK  TUBE 


In  studies  such  as  the  one  presently  under  consideration,  where  the 
thermodynamic  state  of  a  slug  of  high  purity  gas  must  be  accurately  known, 
shock  tubes  have  often  been  successfully  employed  (  References  1-3)  .  The 
shock  tube  offers  several  unique  advantages  over  other  devices  which  produce 
high  temperature  tlows.  Some  of  these  advantages  are:  1)  the  state  of  the  gas 
in  the  shock  tube  can  be  precisely  determined  from  a  knowledge  of  properties 
of  the  gas  and  the  velocity  of  the  shock  as  it  traverses  the  tube;  thus  it  is  not 


lecessary  to  disturb  the  gas  sample  with  probes  in  order  to  deter 


mine  its  state; 


2.)  the  purity  of  the  test  gas  is  controllable,  as  it  never  makes  contact  with  an 
electrode  or  any  high  temperature  surface;  and  3)  as  long  as  the  shock  velocity 
and  conditions  ahead  of  the  shock  are  uniform,  the  gas  sample  produced  will 
also  be  uniform,  since  heat  transfer  to  the  tube  walls  is  negligible  during  the 
short  time  involved  in  any  test. 
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In  the  current  study,  all  the  above-mentioned  advantages  of  the  shock 
tube  are  needed.  Electron  density  and  plasma  frequency  of  air  are  strong 
functions  of  the  thermodynamic  state  and  the  impurity  level  of  the  sample, 
and  both  must  be  accurately  known  and  controlled.  Uniform  gas  samples  are 
needed  tor  the  calibration  of  the  microwave  equipment. 

As  discussed  in  Section  3  of  this  chapter,  it  was  desired  to  construct 
a  shock  tube  to  produce  a  slug  of  hot  air,  or  plasma,  about  18  inches  in  length 
and  -4  inches  in  diameter,  with  equilibrium  electron  densities  n  between  10^ 
and  10  electrons/cc ,  i.e.  plasma  frequencies  v  between  3  x  10?  and  3  x  10^ 
c.p.  s.  lioth  uniform  and  nonunitorm  electron  density  samples  were  required, 
with  the  nonunitorm  samples  having  a  known  axial  electron  density  gradient. 

As  a  preliminary  step  in  determining  the  shock  tube  needed  for  the  study, 
the  equilibrium  electron  density  behind  a  shock  was  calculated  for  various  shock 
Mach  numbers  Mg  in  air  at  room  temperature  and  various  pressures.  Since  it 
was  desired  to  avoid  high  pressures  in  the  test  section  to  facilitate  the  microwave 
detector  installation,  it  was  decided  to  design  the  experiment  for  an  initial  pres¬ 
sure  p,  of  0  .  0  1  atmospheres.  This  pressure  level  is  easily  achieved,  and 
results  in  a  maximum  test  section  pressure  (  after  shuck  reflection)  of  8.  5 
atmospheres.  The  results  of  these  electron  density  calculations,  based  on  data 
from  References  4  and  5,  are  shown  in  Figure  1. 

In  order  to  produce  a  significant  axial  variation  of  electron  density  in  the 
test  gas,  it  was  thought  necessary  to  create  an  axial  temperature  gradient  behind 
the  shock  moving  down  the  tube,  since  electron  density  of  a  gas  is  primarily  a 
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function  of  temperature.  The  most  convenient  method  of  creating  a  temperature 
gradient  behind  the  incident  shock  is  by  varying  the  temperature  in  front  of  the 
shock.  This  will  cause  both  the  shock  velocity  and  shock  Mach  number  to 
change  as  the  wave  passes  through  the  initial  temperature  gradient  and,  by 
careful  selection  of  the  initial  shock  Mach  number,  initial  temperature,  and 
initial  temperature  gradient,  it  should  be  possible  to  produce  the  desired 
gradient  behind  the  shock. 

1  By  means  of  a  simplified  analysis,  the  results  of  various  Mach  number  , 
temperature,  and  temperature  gradient  combinations  have  been  estimated,  and 
the  most  promising  of  these  hav  e  then  been  accurately  calculated  by  means  of  a 
real  gas  characteristics  program  for  the  Bendix  G- 1  5  computer. 

The  simplified  analysis  is  performed  as  follows:  We  choose  values  of 

Mc .  ,  T,  and  AT.  or  T,  (for  meanings  of  symbols,  see  nomenclature  list!  . 

tj-  1  If 

We  assume  that  hydrogen  at  room  temperature  is  the  driver  gas,  and  air  at 
0.0  1  atm.  is  the  driven  gas.  It  is  well  known  that  the  diaphragm  pressure 
ratio  pg/pj  is  a  function  of  M  ,  a(  ,  and  a^.  We  can  thus  compute  P.j/Pj  for 
the  aj  corresponding  to  Tj^,  and  for  Mg..  This  same  pressure  ratio  is  then 
used  to  compute  MS(.  assuming  a^  corresponds  to  Tj(-  We  now  have  two  different 
shock  Mach  numbers  which  correspond  to  the  two  different  initial  temperatures 
at  the  ends  of  the  initial  temperature  distribution.  By  using  References  4  and  5 
we  can  determine  T^.,  T £  .,  p^,  an^  nej  anc'  nep  frum  these.  Although  this 

analysis  does  not  give  details  of  the  electron  density  gradient,  it  does  predict 
the  maximum  variation  of  electron  density  we  can  expect  in  the  plasma  behind 
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the  incident  shock.  When  the  results  of  such  an  analysis  were  satisfactory, 
a  detailed  one -dimensional  unsteady  characteristics  solution  was  obtained  for 
the  same  values  of  Ms^  ,  Tj.  ,  and  Tj^.  The  computer  program  was  a  modified 
version  of  the  one  described  in  Reference  6. 

The  results  of  the  preliminary  analysis  showed  that  a  satisfactory 

gradient  for  a  typical  test  <  se  might  be  established  by  choosing  M s ^  =  5.5, 

T,  =  450  K,  and  T>  =  30  0  K,  withp,  =  0.0  1  atm.  These  numbers  were 
i  1  f  1 

then  fed  intu  the  Hendix  computer  in  a  nondimensional  form  which  specified 
that  the  initial  temperature  in  the  driven  tube  decreased  linearly  with  distance 
along  the  tube.  The  length  over  which  the  temperature  change  occurred  was 
arbitrary  in  the  nondimensional  coordinates  , 

From  the  numerical  solution  of  the  problem  it  was  found  that  the  shock 

slows  down  as  it  traverses  the  region  of  decreasing  initial  temperature,  but  the 

shock  Mach  number  increases  due  to  the  decreasing  value  of  sound  speed  ahead 

of  the  shock.  The  final  shock  Mach  number  was  6,  Z.  It  was  also  found  that 

I  t.  =  .1570  K ,  Ty  .  =  1Z0  0  °K,  r  ,  =  1,5x10  and  n  =  1,1x10  el<;ctrons/cc, 
c\  -i  ef 

For  reasons  that  will  be  discussed  below,  it  was  thought  convenient  to 
make  the  driven  tube  10  feet  long,  with  the  lest  section  extending  this  by  several 
feet.  If  we  heat  the  tube  as  shown  in  Figure  Z,  with  the  initial  temperature 
gradient  coming  between  the  35-foot  and  40-foot  stations  ,  and  if  we  consider 
the  instant  when  the  shock  is  43.5  feet  from  the  diaphragm  station,  i.  e.  5.5 
feet  into  the  test  section,  then  the  plasma  temperature  and  electron  density  in 
the  test  section  will  be  as  shown  in  Figure  3.  It  can  be  seen  that  the  resulting 
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electron  density  gradient  falls  within  the  desired  range  of  values. 

It  should  be  noted  that  the  above  analysis,  while  based  on  real  gas 
computations,  does  not  take  into  account  shock  velocity  nonuniformities  due 
to  such  phenomena  as  boundary  layer  growth  and  finite  diaphragm  opening  time. 
Furthermore,  no  accurate  theory  exists  which  can  predict  these  nonuniformities. 

In  an  actual  test,  we  plan  to  avoid  the  problem  of  accurate  prediction  of  shock 
velocity.  This  will  be  done  by  measuring  the  shock  velocity  at  eight  stations  along 
the  tube.  These  measurements,  along  with  a  knowledge  of  the  initial  temperature 
distribution  in  the  driven  tube,  will  permit  an  accurate  calculation  of  the  actual 
electron  density  gradient  behind  the  shock.  The  computer  will  perform  these 
calculations. 

In  order  to  correctly  calculate  electron  density  using  velocity  measurement?, 
these  measurements  must  be  made  with  great  precision.  Analysis  has  shown  that 
the  error  in  electron  density  is  given  by  the  following  relation: 


d  n. 


d  In  n 


d  T 
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(lus  , 


9  f  n  n. 
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us  J  dP'Po  /T  Po  \  T1 
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where  dx/Xis  defined  as  the  per  cent  error  in  the  assumed  value  of  X.  For  our 
range  of  interest  (  T,  =  110  0  -  3000°K)  this  becomes: 
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J  I 


Thus,  if  we  can  measure  u  with  0  .  1%  accuracy  we  will  still  have  an  error  of  more 


than  5%  in  electron  density.  The  error  of  a  few  degrees  that  we  may  make  in  T 
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will  have  little  effect  on  the  error  in  n^.  We  plan  to  accomplish  the  velocity 

measurement  by  using  eight  pair  of  ionization  probes  distributed  along  the 

driven  tube  and  test  section.  Each  pair  of  probes  will  constitute  one  velocity 

measuring  station,  with  the  two  probes  about  two  feet  apart.  The  ionization 

behind  the  shock  wave  will  trigger  each  probe  in  turn,  and  the  probes  will  act 

as  shock  time  -  ot'-ar rival  indicators.  Special  circuits  have  been  designed  to 

-  6 

measure  and  record  the  outputs  with  the  required  accuracy  (0.2  x  10  seconds) 

The  primary  factor  which  was  used  to  fix  the  length  of  the  shock  tube  was 

the  length  of  the  plasma  slug  required.  If  ideal  gas  relations  are  used,  the 

u,  1 

plasma  length  per  unit  length  of  driven  tube  is  given  by  1  -  .  This 

will  be  reduced  due  to  real  gas  effects  and  attenuation.  For  5.5  Mg  6 
we  estimate  real  plasma  lengths  of  about  0,0  8  feet  per  foot  of  driven  tube. 

Now,  for  our  test  we  will  require  about  3.5  feet  of  plasma,  and  thus  we  have 
chosen  to  use  a  40  -  foot  driven  tube  plus  a  test  section  about  4  feet  long.  The 
driver  is  10  feet  long.  The  entire  shock  tube,  including  driver  section,  has  a 


nominal  inside  diameter  of  three  inches. 


The  temperature  distribution  is  obtained  by  winding  insulated  nichrome 
wire  helically  around  the  driven  tube.  The  wire  is  heated  by  passing  current 
through  it,  and  is  divided  into  15  individually  controlled  sections  which  allow 
the  establishment  of  any  desired  initial  temperature  distribution.  The  entire 
driven  tube  is  thermally  insulated  from  the  surrounding  air  and  supporting  stand. 

Heat  from  the  tube  wall  will  be  conducted  to  the  low-pressure  air  in  the 
driven  tube,  and  the  initial  temperature  distribution  in  the  air  should  follow  closely 
that  of  the  tube  wall.  Convection  is  not  expected  to  play  a  major  role  in  the  transfer 


process. 
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Although  it  was  originally  thought  that  our  tests  would  be  performed 
in  the  GASL  shock  tunnel,  we  soon  realized  that,  due  to  the  moderate  values  of 
temperature  and  pressure  required,  it  would  be  more  economical  and  convenient 
to  construct  a  new  shock  tube  with  relatively  thin  walls.  Much  less  power  is 
required  to  heat  this  new  tube,  which  is  made  of  stainless  steel  pipe,  to  the 
desired  initial  temperature  distribution,  and  the  cost  of  building  the  tube  is 
considerably  less  than  that  of  extending  the  shock  tunnel  tube  to  a  length  of  40  feet. 

In  the  body  of  shock  tube  literature  available  no  reference  is  to  be  found 
of  any  facility  using  this  type  of  differential  preheating  of  the  driven  gas.  Thus 
we  believe  we  have  developed  a  unique  and  versatile  research  tool. 

From  the  values  of  shock  Mach  number  and  initial  pressure  to  be  used, 
we  can  calculate  the  hydrogen  pressure  needed  in  the  driver.  This  is  of  the 
order  of  10  atmospheres;  a  nonmetallic  diaphragm  is  used  to  contain  this.  In 
order  to  control  the  shock  Mach  number  produced  we  will  charge  the  driver  and 
driven  tubes  to  the  desired  pressures  using  a  diaphragm  which  will  not  burst  at 
these  pressures.  Then  a  wire  in  contact  with  the  diaphragm  will  be  heated 
electrically,  and  thus  weaken  the  diaphragm  material  until  it  bursts.  The  amount 
of  heating  required  to  produce  a  good  burst  should  be  small.  The  diaphragm  is 
thermally  insulated  from  the  heated  driven  tube  by  phenolic  rings. 

In  an  attempt  to  estimate  the  effects  of  boundary  layer  on  the  flow  in  the 
shock  tube,  the  distance  behind  the  shock  at  which  transition  occurs  has  been 
calculated.  Assuming  a  transition  Reynolds  number  of  6  x  10  ^  based  on  the 
distance  a  particle  has  moved  after  the  shock  has  passed  over  it,  we  found  that 
transition  occurs  24  inches  behind  the  shock  in  the  test  section.  Thus  we  will  have 
about  two  feet  of  laminar  boundary  layer. 


57 


The  electron  density  profile  in  the  laminar  boundary  layer  ha 


s  been 


computed,  using  data  from  Refe 


rence  7  with  the  equilibrium  value  for  plasma 


in  Figure  4.  We  sec 


temperature.  The  results  of  this  calculation  ar^shown  i 
that  a  Mach  6  shock,  moving  into  air  at  300°K  and  0 . 0  1  atmospheres  (approxi¬ 
mately  the  conditions  existing  in  the  test  section)  ,  produces  a  boundary  layer  in 
which,  two  feet  behind  the  shock,  the  electron  density  has  decreased  10%  from 
tree  stream  at  0. 18  inches  from  the  tube  wall.  We  can  define  a  core  area, 
within  which  the  electron  density  is  no  more  than  10%  from  its  free  stream  value, 
and  the  ratio  of  core  area  to  total  area  will  vary  as  in  Figure  5. 

In  the  course  of  our  tests,  an  antenna  at  the  end  of  the  test  section  will 
be  radiating  electromagnetic  eaves  upstream,  through  the  oncoming  shock,  and 
into  the  hot  plasma  behind  the  shock.  To  prevent  this  radiation  from  continuing 
upstream,  reflecting  off  the  end  of  the  driver,  and  coining  back  downstream 
interfere  with  the  tests,  we  have  designed  an  electromagnetic  trap  which,  when 
installed  in  the  tube,  will  transform  the  test  section  into  an  electromagnetic 
cavity.  This  trap  will  consist  of  one  or  more  metallic  ribbons  stretched  across 
the  tube  between  the  test  section  and  the  driven  tube.  The  ribbons  will  be  0.00  1 
inches  thick  with  a  chord  of  about  0.25  inches,  and  will  be  in  tension  to  prevent 


fluttering.  The  flow  in  the  test  section  i 


s  not  expected  to  be  affected,  since  the 


ribbons  will  be  several  inches  upstream  of  the  re  it 


gion  of  interest. 
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3.  EXPERIMENTAL  INVESTIGATION  OF  THE  PLASMA  PROPERTIES 


The  preceding  section  discussed  the  aerodynamic  requirements  and 
properties  of  the  GASL  shock  tube  built  for  this  experiment.  In  this  section, 
we  discuss  the  experimental  facility  to  deti  rmine  the  electromagnetic  properties 
of  a  nonuniform  plasma.  The  use  of  a  shock  tube  to  generate  a  plasma  is  highly 
advantageous  for  the  reasons  previously  mentioned.  However,  using  the  plasma 
region  behind  the  shock  wave  imposes  severe  requirements  on  the  experimental 
technique  to  measure  the  electromagnetic  properties  of  the  ionized  slug.  The 
primary  difficulty  is  due  to  the  short  duration  that  the  plasma  region  is  properly 
located  in  the  electromagnetic  test  section.  For  a  uniform  plasma  such  as  that 
generated  in  a  conventional  shock  tube,  tin  running  times  are  on  the  order  of 
100  p.sec.  However,  for  an  axially  nonuniform  plasma,  it  is  clear  that  the  electro¬ 
magnetic  properties  are  constantly  changing  while  the  slug  is  passing  through 
the  test  section.  A  number  of  interesting  phenomena  may  occur  as  the  slug 
enters  the  testing  region.  Such  phenomena  may  have  very  short  transients. 

For  these  reasons,  it  was  deemed  necessary  that  the  detection  systt-m  have  a  high 
time  resolution  ability. 

The  second  condition  to  be  satisfied  by  the  ionized  gas  is  given  by  the 
maximum  value  of  the  electron  density  at  which  the  e lectromagnetic  wave  can 
propagate  through  the  plasma.  A  suitable  compromise  between  mechanical 
feasibility  of  the  shock  tube  and  a  convenient  frequency  of  the  electromagnetic 
wave  suggested  a  tube  diameter  of  3,94  inches.  This  diameter  is  maintained 
for  both  the  driven  and  test  sections.  In  the  cylindrical  wave  guide  corresponding 
to  the  test  section,  the  T  En  mode  was  selected  for  the  electromagnetic  wave. 
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This  mode  allows  the  use  of  a  minimum  frequency  in  the  range  of  3x10  cycles  per  sec. 
On  the  other  hand,  this  frequency  is  lower  than  the  cut-off  frequency  of  the 
T  Mq1  propagation  mode.  The  frequency  of  the  electromagnetic  wave  determines 
the  maximum  value  of  the  electron  density  in  the  ionized  gas  which  corresponds 
to  approximately  10^  electrons  per  cc. 

The  minimum  electron  concentration  is  determined  by  the  minimum  detectable 
signal  from  the  electromagnetic  probes.  With  the  present  technique,  it  is 

7 

projbable  that  the  minimum  detectable  electron  density  is  approximately  10 
electrons/cc. 

Since  it  is  the  purpose  of  this  investigation  to  study  the  properties  of  non- 
uniform  plasma,  the  electron  density  in  the  extreme  case  should  be  rapidly 
changing  in  a  distance  comparable  to  one  wavelength.  By  establishing  a  tem¬ 
perature  gradient  in  the  tube  as  discussed  in  the  previous  section,  an  electron 
gradient  as  strong  as  one  order  of  magnitude  every  two  wavelengths  can  be  realized. 

Another  difficulty  in  using  the  plasma  generated  behind  the  shock  wave  is  that 
the  electron  density  is  not  radially  uniform.  This  lack  of  uniformity  leads  to 
an  average  electron  concentration  which  is  lower  than  the  electron  density  at 
the  axis  of  the  tube.  Since  the  electrodynamic  properties  of  the  plasma  can  only 
be  conveniently  measured  by  considering  the  whole  tube  as  a  waveguide  uniformly 
filled  with  a  dielectric  medium,  the  radial  plasma  nonuniformity  must  be 
minimized.  Since  the  radial  plasma  nonuniformity  is  closely  related  to  the 
aerodynamic  boundary  layer,  great  care  was  taken  in  the  design  and  construction 
of  all  parts  of  this  facility  to  minimize  the  boundary  layer  growth  as  the  shock 
wave  proceeds  down  the  tube.  The  effect  of  the  radial  nonuniform  plasma 
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distribution  is  being  analyzed  from  a  theoretical  point  of  view  as  contained 
in  Chapter  III.  The  aforementioned  conditions  determine  the  basic  design 
considerations  for  this  facility.  The  major  elements  of  this  design  are  discussed 
be  low. 

In  Figure  7,  the  end  of  the  driven  part  of  the  shock  tube  is  shown.  A 
microwave  cavity  is  formed  between  the  ribbon  section  and  the  tuning  slug  in 
the  antenna  section.  The  microwave  power  in  the  range  of  3  kilomegacyeles 
is  fed  to  the  antenna  which  is  stretched  across  the  tube  approximately  l/4 
wavelength  from  the  tuning  slug.  The  complete  microwave  system  is  outlined 
in  block  form  in  Figure  6. 

The  detection  assembly  is  located  three  wavelengths  from  the  antenna 

section  to  avoid  the  effects  of  the  attached  field  in  the  antenna  region.  Ten 

I 

detectors  are  mounted  at  quarter  wavelength  intervals  for  the  median  waveguide 
wavelength.  Each  detector  is  a  modified  crystal  with  an  individually-adjustable 
capacitive  coupling.  The  probe  is  mounted  flush  with  the  wall  of  the  tube  to 
minimize  boundary  layer  growth.  The  coupling  between  a  probe  and  the  wave¬ 
guide  must  be  kept  at  a  minimum  value  to  avoid  a  large  perturbation  of  the  field 
distribution  in  the  waveguide. 

In  order  to  record  data  from  the  detectors,  special  ultra  high  gain  amplifiers 
are  needed.  Since  it  is  desirable  to  observe  the  short  transient  effects  as  the 
shock  passes,  these  amplifiers  were  required  to  have  a  band  pass  of  7  megacycles. 
Since  the  data  of  each  detector  must  be  recorded  independently,  it  was  deemed 
best  that  individual  oscilloscopes  and  cameras  be  employed  for  each  detector. 

Since  no  commercially  available  scope  with  the  desired  band  pass  and  sensitivity 
was  available  within  reasonable  cost,  the  ten  scope  units  were  designed  and 


61 


built  by  the  GASL  facility.  The  instrument  design  incorporates  a  precision 
attenuator  and  a  variable  bandpass  switch  to  accurately  determine  the  short  and 
medium  drift  in  the  microwave  generator.  All  of  the  scopes  are  driven  with 
a  common  sweep  and  trigger  source.  The  specification  on  each  scope  is 
as  follows: 

Deflection  sensitivity-from  100  pv/cm  to  l  v/cm  in  calibrated  steps. 

Rise  time  -  .050  p  sec. 

Bandpass  -  l /Z  to  7  me  adjustable. 

Noise  -  input  open  circuited  -  ZO  pv  RMS. 

-  input  short  circuited  -  10  pv  RMS. 

It  should  be  noted  that  each  detector  is  AC  coupled  to  the  amplifiers  so  that 
only  the  change  in  microwave  signal  is  displayed  on  the  scope. 

In  Figu  re  8  a  field  distribution  is  shown  typical  of  that  encountered  if  the 
plasma  has  an  axial  gradient.  Note  that  the  wavelength  has  shifted  and  that  the 
signal  is  partially  attenuated. 

For  comparison,  the  field  distribution  for  the  waveguide  without  plasma  is 
superimposed.  The  recorded  data  during  the  plasma  flow  will  determine  the 
instantaneous  departure  of  the  field  distribution  with  respect  to  the  stationary 
distribution  prior  to  each  experiment.  This  data  is  recorded  on  a  continuous 
time  basis  during  the  running  time  of  the  tube.  Hence,  for  each  detector,  a 
history  of  the  difference  between  the  zero  plasma  case  and  the  incident  plasma 
is  made.  By  comparing  the  ten  histories  thus  recorded,  the  field  distribution 
and  hence  the  wavelength  and  attenuation  can  be  determined  for  all  times  during 
the  test.  A  particularly  interesting  situation  may  be  found  when  the  electron 
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density  increases  behind  the  shock  reaching  the  critical  value  which 
corresponds  to  a  plasma  frequency  equal  to  the  angular  frequency  of  tne 
electromagnetic  wave.  The  position  in  the  ionized  layer  of  this  critical 
density  corresponds  to  the  transition  from  the  propagation  region  to  a  region 
of  complete  attenuation  of  the  electromagnetic  wave.  The  passage  of  this 
transition  plane  by  the  position  of  a  detector  leads  to  a  drastic  change  in  the 
time -dependent  signal  obtained  by  the  detector.  We  point  out  that  it  is  of 
major  interest  to  analyze  the  correlation  between  the  theoretical  calculation 
and  the  experimental  results  across  this  transition  region.  This  is  due  to  the 
fact  that  the  position  of  this  transition  region  in  the  plasma  may  be  strongly 
affected  by  even  a  small  departure  of  the  actual  electromagnetic  properties 


of  the  plasma  from  the  theoretical  model. 

At  the  present  date,  all  of  the  mechanical  work  on  the  shock  tube 
has  been  completed.  In  addition,  all  of  the  microwave  equipment  and  the  test 
section  has  been  finished.  The' data  recording  system  composed  ol  the  detec¬ 
tors,  scope  amplifiers,  scopes  and  sweep  source  are  complete  except  for 
calibration.  In  like  manner,  the  heating  and  control  system  for  establishing 
a  variable  gradient  in  the  tube  is  finished.  The  only  remaining  pretest  work  on  the 
facility  is  principally  a  careful  calibration  of  all  components.  It  is  expected 
that  preliminary  test  runs  will  be  started  shortly. 
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SYMBOLS  -  CHAPTER  V 


a  speed  of  sound 

A  tube  cross  section  area 

M  shock  Mach  number  :  u  /a, 

s  s 

ne  electron  density 

p  pressure 

T  temperature 

AT1  Tj.  -  Tlf 

u  velocity 

X  any  variable 

x  axial  distance,  either  from  diaphragm  or  behind  shock 

y  normal  distance  from  wall 

Vp  plasma  frequency 

p  density 

Subscripts 

"  ( T  =  275.16  K 

/  o 

o  reference  conditions  | 

Pj  =  1  atm,  p  o  =  i . 288  x  1  0  grn/ cc 

1  initial  conditions  in  driven  tube 

2  conditions  in  plasma  behind  shock  wave 

4  initial  conditions  in  driver  tube 

i  conditions  at  diaphragm  end  of  region  of  interest 

f  conditions  at  test  section  end  of  region  of  interest 

s  pertaining  to  shock  wave 


65 
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FIGURE  6.  MICROWAVE  SYSTEM 
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FIGURE  7.  LONGITUDINAL  CROSS  SECTION  OF  TEST  CHAMBER 
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